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Nomenclature 


A  =  spherical  shell  surface  area  element  in  the  plume 

Ei  =  voltage  exchange  parameter 

E2  =  mass  exchange  parameter 

F  =  thrust  density  vector  in  the  plume 

( j  =  normalized  ion  mass  flow  fraction  of  /th  ion  species 

&  -  Faraday  constant,  96,485  coulombs/mol  of  charge 

/(v)  =  velocity  distribution  function  of  ions  and  neutrals 

g  =  Earth’s  gravitational  constant  at  sea  level,  9.806  m/s2 

Id  =  anode  discharge  current 

I-Beam  =  integrated  beam  current 

I  Axial  =  axial  component  of  beam  current  parallel  to  thruster  centerline 

Isp  =  specific  impulse 

j  =  propellant  charge  state  index,  0,  1,  2,  3,  etc.  for  Xe°,Xe+1,Xe+2,Xe+3 

J(0)  =  current  density  in  the  plume  at  angular  position  9 

rh  =  anode  propellant  mass  flow  rate,  where  X  nr,  =  riifor  j= 0  to  the y'th  ion  species 

m  (9)  =  anode  propellant  mass  flow  rate  at  angular  position  9 

mi  =  mass  flow  rate  of  ions,  where  £  m;  =  mi  for  j=  1  to  the y'th  ion  species 

m;  =  mass  flow  rate  of y'th  propellant  species 

M  =  molecular  weight  of  propellant  (Xe  =  0. 13 13  kg/mol) 

Pd  =  discharge  power  to  anode  =  Vd  Id 

Pjet  =  jet  power 

Pmin  =  minimum  power  loss  to  sustain  ionization 

Q  =  average  charge  of  propellant  ions 

q  =  unit  of  charge,  1.609  x  10'19  Coulombs 

r  =  fraction  of  electron  current  to  the  anode,  electron  recycle  fraction 

(1-r)  =  current  utilization,  fraction  of  ion  beam  current  =  IBeam  /  Id 

rmin  =  minimum  fraction  of  electron  current  required  to  sustain  discharge 

T  =  component  of  thrust  vector  directed  along  thruster  centerline 

Va  =  most  probable  ion  acceleration  voltage 

Vd  =  anode  discharge  voltage 

Wj  =  exit  speed  of  y'th  species 

V  =  average  exit  velocity  of  the  VDF  over  velocity  space  dv  at  angular  position  9 

V  (0)  =  radial  component  of  V  in  hemispherical  coordinates  from  thruster  centerline  at  angular  position  9 

v,  v2  =  average  propellant  velocity,  squared  propellant  velocity 

vT,  vj2  =  average  ion  velocity,  squared  ion  velocity 

y j  -  normalized  speed  ratio  of  the  y'th  species  =  vy  /lv ,  I 

Z j  =  ion  charge  state  =  1,  2,  3  for  Xe+1,  Xe+2,  Xe+3 

p  =  fractional  loss  of  acceleration  potential 

( 1  -p)  =  voltage  utilization,  Va  /  Vd 

SVj  =  acceleration  potential  of y'th  ion  species 

e  j  =  ionization  potential  of  /th  ion  species  (12,  33,  65  eV  from  neutral  ground  state  for  Xe) 
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r|E  =  energy  efficiency  =  (l-(3)  (1-r) 

r|A  =  anode  thrust  efficiency 

9  =  angular  position  in  the  plume,  9=90°  on  thrust  axis 

X  =  plume  momentum  divergence  half-angle,  X=0°  on  thrust  axis 

p  =  radial  coordinate  in  the  spherical  coordinate  system 

<f>m  =  mass  utilization,  ion  mass  flow  fraction  at  exit,  Y,  m;  /  rnfor  j=  1  to  the  jth  ion  species 

<f>n_g  =  neutral  gain  utilization  resulting  from  the  bulk  and  thermal  speed  of  neutral  propellant 

<J>P  =  propellant  efficiency,  losses  due  to  dispersion  of  the  VDF  =  <t>m  <f>N_G  3>q 

<f>q  =  charge  utilization 

X  =  output  moles  of  ion  charge  per  input  moles  of  propellant  =  <I>m  Q  =  IBeam  M l(  rii  &) 

*PB  =  beam  efficiency,  divergence  loss  component  of  thrust 

£1,  =  normalized  ion  species  current  fraction  =  f ;  Z,  /  Q 

Q  =  normalized  ion  species  number  density  fraction  ~  f*,  Z,"1/2 

<  >,„  =  mass  weighted  average  quantity  in  the  plume  (0<9<ji) 

<  >mv  =  momentum  weighted  average  quantity  in  the  plume  (0<9<ji) 

<  >j  =  charge  flux  weighted  average  quantity  in  the  plume  (0<9<7i) 


I.  Introduction 

Standardization  of  experimental  methods,  facilities,  diagnostic  apparatus,  and  efficiency  analysis  has 
previously  been  proposed  in  the  electric  propulsion  community.1  An  efficiency  architecture  is  presented  that  is 
based  on  a  consistent  set  of  definitions,  conservation  of  energy,  conservation  of  mass  and  charge,  and  Newton’s 
Second  Law.  The  relationships  between  Hall  thruster  performance  data,  telemetry  (Va,  Ia,  and  m),  and  plume 
measurements  are  summarized  and  new  correlations  are  developed. 

In  this  paper,  we  present  the  basis  for  analytical  separation  of  anode  thrust  efficiency  into  the  product  of 
energy  efficiency,  propellant  efficiency,  and  beam  efficiency,  which  are  less  than  unity  under  all  operating 
conditions.  Separating  utilization  efficiencies  in  this  manner  isolates  performance  losses  in  terms  of  energy  losses 
that  lead  to  Joule  heating,  dispersion  of  the  jet  velocity  distribution  function  (VDF),  and  loss  of  thrust  due  to  plume 
divergence.  While  these  processes  are  physically  coupled,  the  effects  on  performance  may  be  mathematically 
isolated. 

Historical  perspective  of  early  ion  engine  analyses,  methodologies  developed  in  the  former  Soviet  Union, 
and  a  brief  description  of  contemporary  efficiency  models  are  discussed  in  relation  to  the  presented  architecture. 
The  architecture  is  distinctive  in  formulating  the  axial  component  of  thrust  as  the  product  of  mass-weighted  and 
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momentum-weighted  quantities,  and  introduces  a  new  definition  of  propellant  utilization  that  incorporates  multiple 
ion  species  and  the  effects  of  neutral  propellant.  Thruster  performance  parameters,  including  thrust  to  power  (T/P) 
and  specific  impulse  (Isp),  are  formulated  in  terms  of  utilization  efficiencies  to  provide  a  higher  level  of  confidence 
in  comparisons  between  plume  measurements  and  thrust  measurements.  Experimental  thrust  and  plume 
measurements  from  a  6  kW  laboratory  model  Hall  thruster  are  studied  to  illustrate  the  fidelity  of  the  proposed 
analytical  model  and  to  demonstrate  the  utility  of  combining  results  from  a  set  of  diagnostics  in  determining  Hall 
thruster  performance.  Power  losses  are  studied  using  energy  efficiency  and  compared  to  the  minimum  power 
required  to  sustain  the  discharge  ion  species  composition. 

II.  Historical  Perspective  and  Recent  Efforts 

Investigations  of  ion  acceleration  and  electron  transport  using  Hall  thruster  technology  began  in  the  early 
1960s  in  the  United  States  and  former  Soviet  Union.2,3'4'5,6’7'8'9'10'11'12  The  focus  of  US  electric  propulsion  research 
shifted  primarily  to  ion  engine  technology  in  the  early  1970s13,  whereas  investigations  in  the  USSR  continued  Hall 
thruster  advancements  throughout  the  following  decades.141516 

Analytical  factorization  of  ion  engine  thruster  efficiency  was  described  as  early  as  1975  by  Masek  et  al.  in 
a  review  of  ion  engine  performance.17  Thrust  efficiency  was  factored  into  the  product  of  energy  efficiency  and 
propellant  utilization  efficiency  using  Newton’s  Second  Law,  with  terms  accounting  for  losses  due  to  doubly 
charged  ions  and  beam  divergence.  However,  the  derivation  was  condensed  and  an  explanation  for  the  treatment  of 
beam  divergence  and  multiply  charged  ions  was  not  presented.  The  methodology  was  not  widely  adopted  or  cited  in 
the  ion  engine  community  and  is  absent  in  modem  analysis. 

In  the  1990s,  the  manifestation  of  Russian  Hall  thruster  technology  in  the  Western  and  Japanese  spacecraft 
communities  catalyzed  a  resurgence  of  Hall  thruster  research.  '  '  ~  ’  Technology  transfer  that  began  in  the 
early  1990s  brought  invaluable  benefits  and  advancements  to  Western  Hall  thruster  development,  but  much  of  the 
earlier  Soviet  progress  in  Hall  thmster  research  that  was  published  in  the  Russian  language  was  not  translated.  The 
limited  availability  of  translated  documents  from  this  extensive  literature  inevitably  impeded  the  transmission  of 
knowledge  and  progress.  As  a  result,  several  important  contributions  from  Russian  research  have  not  been  widely 
disseminated  in  the  West,  including  the  analytical  factorization  of  anode  thrust  efficiency. 
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Factorization  of  Hall  thruster  efficiency  was  outlined  in  a  manual  on  stationary  plasma  engines  by  Belan, 
Kim,  Oranskiy,  and  Tikhononv  from  the  Kharkov  Aviation  Institute  in  1989. 24  This  seminal  document  was  later 
referenced  and  the  performance  methodology  summarized  by  Bugrova  et  al.  without  explicitly  stating  the  nature  of 
the  earlier  derivation. 2:1  Kim’s  highly  cited  1998  paper  on  processes  that  determine  Hall  thruster  efficiency26  and 
other  contemporary  publications  ’  '  continued  the  elaboration  of  the  Kharkov  Aviation  Institute  methodology. 
These  approaches  showed  similarities  to  the  analysis  presented  here,  and  demonstrated  relationships  between 
experimental  variables  and  plasma  phenomena  that  affect  thruster  performance. 

A  textbook  by  Grishin  and  Leskov  analyzed  thruster  performance  based  on  energy  efficiency  and  thrust 
efficiency.30  The  ratio  of  thrust  efficiency  to  energy  efficiency  was  asserted  to  describe  the  velocity  dispersion  in 
magnitude  and  direction.  This  ratio  will  be  derived  from  first  principles  in  the  following  section,  where  it  will 
become  apparent  that  the  ratio  is  equivalent  to  beam  efficiency  and  propellant  efficiency,  which  capture  losses 
resulting  from  plume  divergence,  incomplete  ionization,  and  the  production  of  multiply  charged  ions. 

The  term  propellant  utilization  efficiency  has  been  used  in  at  least  three  different  ways  in  the  literature. 
Most  commonly  as  (1)  the  ionization  fraction,  as  (2)  the  fraction  of  momentum  carried  by  ions14,  or  as  (3)  the  ratio 
of  output  moles  of  charge  to  the  input  moles  of  propellant31,  which  is  equivalent  to  the  product  of  the  ionization 
fraction  and  average  ion  charge.  While  the  definitions  of  voltage  and  current  utilization  are  standardized,  a 
consistent  description  of  propellant  utilization  has  not  emerged. 

In  the  post- 1990  era,  numerous  efficiency  analysis  frameworks  and  modifications  have  been  proposed, 32'33, 
34,35,36,37,38,39,40  jncju(jjng  notable  studies  by  Komurasaki41,  Ahedo42,  and  Hofer43'44.  The  difference  between  the 
proposed  architecture  with  past  ion  engine  studies,  the  analytical  factorization  in  Russian  literature,  and  the  post- 
1990  methodologies  varies  in  each  case.  Dissimilarities  arise  due  to  varying  levels  of  completeness  regarding  the 
treatment  of  multiply  charged  ions,  beam  divergence,  and  the  effect  of  neutral  propellant  on  dispersion  of  the  jet 
VDF.  In  some  instances,  the  effects  of  these  loss  mechanisms  are  neglected  or  inserted  as  a  utilization  efficiency 
without  rigorously  factoring  the  term  from  the  anode  thrust  efficiency  equation.  Utilization  efficiencies  in  the 
performance  model  presented  in  this  paper  are  analytically  separated  from  the  thrust  efficiency  and  are  formulated  to 
minimize  the  introduction  of  new  terminology.  The  primary  differences  with  all  previous  methodologies  is 
quantifying  the  effect  of  neutrals  on  dispersion  of  the  VDF  in  the  neutral-gain  utilization  and  decomposing  total 


5 


thrust  into  the  product  of  mass-weighted  and  momentum-weighted  quantities.  Three  appendices  are  included  to 
define  plume  averaged  terms  and  detail  the  derivation  of  propellant  efficiency  and  beam  efficiency. 

This  analytical  tool  is  not  meant  to  predict  Hall  thruster  physics  or  plasma  properties  in  the  same  manner  as 
a  computational  model.  Numerical  simulations  typically  determine  plasma  properties  throughout  the  plume  and/or 
performance  attributes  based  on  a  particle  source  mode.  In  contrast,  the  analytical  model  employs  performance 
measurements  and  bulk  plume  characteristics  to  experimentally  determine  physical  processes  and  relationships  that 
are  difficult  to  measure  directly. 

III.  Hall  Thruster  Performance  Architecture 

A.  Decoupling  Energy  Efficiency,  Beam  Efficiency,  and  Propellant  Efficiency 

The  total  thrust  generated  by  a  Hall  thruster  is  a  primarily  a  function  of  electromagnetic  forces,  and  to  a 
much  lesser  degree  gasdynamic  forces.  Downstream  of  the  primary  ionization  and  acceleration  regions,  the  total 
resultant  thrust  may  be  found  from  Eq.  (1)  through  integration  of  the  thrust  density  vector  throughout  the  plume. 

T/3  =  JJf  •/>  dA  (1) 

This  idealized  description  suffers  from  numerous  experimental  difficulties  and  uncertainties.  Although  Eq. 
(1)  may  be  valid  in  the  space  environment,  facility  effects  inherent  in  ground  tests  result  in  significant  scattering  in 
the  plume  and  potentially  increased  thrust  from  neutral  ingestion.  Measurements  of  thrust  density  throughout  the 
HET  plume  are  extremely  difficult,  but  may  be  viable  using  an  impact  target  plate.45 

The  conventional  inverted  pendulum  thrust  stand  measures  the  component  of  thrust  directed  along  the 
thruster  centerline  axis.  Careful  thruster  alignment  and  uniform  propellant  injection  generally  result  in  an 
axisymmetric  plume  with  negligible  deviations  between  the  resultant  thrust  vector  direction  and  the  thruster 
centerline  axis.  For  an  axisymmetric  plume,  the  steady-state  scalar  component  of  thrust  directed  along  thruster 
centerline  is  formulated  in  Eq.  (2).  The  total  thrust  may  be  factored  into  the  product  of  mass  flow,  mass-weighted 
average  velocity,  and  momentum-weighted  average  divergence  using  the  definitions  listed  in  Appendix  A.  Thrust  is 
often  defined  as  rh  <C>  <cos(X)>  without  mathematically  separating  terms  from  the  thrust  integral. 

71/2 

T  =  27iR2  | rii(9)  v (9)cos(9)sin(9)d9  =  rh  <  v  >m<cos(9)>mv  (2) 

0 
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In  Eq.  (3),  the  standard  definition  of  anode  thrust  efficiency  is  expressed  in  terms  of  thruster  telemetry  and 


measured  thrust.  Anode  efficiency  is  decomposed  into  the  product  of  energy  efficiency,  propellant  efficiency,  and 
beam  efficiency  using  the  definition  of  jet  kinetic  energy,  Newton  s  Second  Law,  and  the  mass-weighted  average 
squared  velocity  as  defined  in  Eq.  (A-5)  of  Appendix  A. 
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Energy  efficiency  in  Eq.  (4)  characterizes  the  conversion  of  input  anode  electrical  energy  to  jet  kinetic 
energy,  and  contains  all  information  about  losses  that  ultimately  appear  as  Joule  heating  of  the  channel  walls  and 
thruster  body,  radiation  from  the  jet,  and  frozen  ionization  losses  in  the  plume.  In  this  analysis,  the  jet  power 
accounts  for  the  energy  of  neutral  propellant,  whereas  most  definitions  include  only  beam  ion  energy. 
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Jet  momentum  losses  due  to  beam  divergence  are  quantified  using  the  momentum-weighted  average  cos(9) 
in  Eq.  (5).  This  formulation  is  similar  to  the  focusing  efficiency  described  by  Kim,  but  is  naturally  expressed  as  a 
momentum-weighted  average  quantity  from  the  formulation  of  thrust. 

'Pb  =<  cos  (0)  >mv  (5) 

Propellant  efficiency  in  Eq.  (6)  is  the  mathematical  relationship  between  the  particle  momentum  and  the  jet 
kinetic  energy.  It  contains  all  loss  information  associated  with  dispersion  of  the  jet  VDF  due  to  incomplete 
ionization  and  the  presence  of  multiple  ion  species  with  widely  varying  velocities.  This  ratio  is  unity  for  100% 
ionization  to  a  single  ion  species. 
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The  factors  affecting  propellant  efficiency,  energy  efficiency,  and  beam  efficiency  will  be  discussed  in  the 
following  sections.  Propellant  efficiency  will  be  decomposed  into  mass  and  charge  utilization  efficiencies,  and  the 
neutral-gain  utilization  is  introduced  to  characterize  the  effect  of  neutral  propellant  on  dispersion  of  the  VDF.  Next, 
energy  efficiency  is  separated  into  the  product  of  voltage  utilization  and  current  utilization  using  %.  Beam  efficiency 
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will  then  be  analyzed  in  terms  of  the  distinction  between  momentum-weighted  divergence  and  current-weighted 
plume  divergence. 


B.  Decoupling  Mass,  Charge,  and  Neutral-Gain  Utilization  from  Propellant  Efficiency 

The  influence  of  non-uniform  velocity  distribution  on  propellant  efficiency  due  to  multiple  ion  species  and 
unionized  propellant  is  expressed  in  Eq.  (7).  Propellant  efficiency  is  separated  into  the  product  of  mass  utilization, 
charge  utilization,  and  neutral-gain  utilization  using  the  definitions  in  Appendix  B.  The  propellant  efficiency  is 
separated  in  this  way  to  maintain  a  consistent  definition  of  mass  utilization  and  charge  utilization  with  previous 
methodologies43,  and  isolate  effects  resulting  from  the  bulk  and  thermal  speed  of  neutrals. 
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The  standard  definition  of  charge  utilization  is  formulated  from  the  first  term  in  brackets  in  Eq.  (7).  Charge 
utilization  is  expressed  in  Eq.  (8)  using  species  current  fractions  as  obtained  from  an  ExB  probe  and  is  unity  for 
ionization  to  a  single  ion  species.  Charge  utilization  is  also  written  using  ion  species  mass  flow  fractions  in  Eq.  (8). 
Ion  mass  flow  fractions  are  defined  according  to  Eq.  (9). 
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It  is  important  to  note  the  distinction  between  ion  current  fractions,  ion  flow  fractions,  and  ion  species 
fractions.  Ion  current  fraction  is  determined  directly  from  ExB  probe  traces.  Ion  species  fractions  are  based  on 
particle  number  density,  and  are  commonly  reported  for  incorporation  in  numerical  simulations.  The  ion  flow 
fraction  is  a  more  suitable  figure  of  merit  for  experimental  plume  studies,  since  it  is  related  to  the  ion  species 
momentum  in  the  plume.  The  average  ion  charge  Q  in  Eq.  (10)  is  expressed  in  terms  of  ion  current  fractions,  ion 


flow  fractions,  and  ion  species  fractions. 
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The  mass  utilization  and  neutral-gain  utilization  are  developed  from  the  second  term  in  brackets  in  Eq.  (7). 
Mass  utilization  is  defined  as  the  ratio  of  ion  mass  flow  rate  to  anode  mass  flow  rate  in  Eq.  (11),  and  is  calculated 
using  beam  current  from  a  Faraday  probe  in  conjunction  with  the  average  ion  charge  from  an  ExB  probe. 
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The  effects  of  neutral  propellant  on  dispersion  of  the  jet  VDF  are  quantified  in  Eq.  (12).  This  term  is 
always  greater  than  one,  and  is  indicative  of  the  thmst  and  energy  gained  due  to  the  speed  of  the  neutrals.  The 
neutral-gain  utilization  is  a  strong  function  of  the  normalized  neutral  speed  and  mass  utilization.  The  normalized 
neutral  speed  is  typically  less  than  0.05  for  most  Hall  thrusters  and  operating  conditions,  which  results  in  a  neutral- 
grain  utilization  of  less  than  1.02.  Neutral-gain  utilization  is  expected  to  increase  for  low  voltage,  high  power 
thruster  operation  and  anode  propellant  injection  with  a  large  bulk  axial  velocity  component.  Minimizing  this 
utilization  is  critical  for  high  performance  operation,  since  it  increases  with  decreased  ionization  and  decreased 
neutral  residence  time  in  the  discharge  channel. 
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The  factorization  of  propellant  efficiency  is  detailed  in  Appendix  B,  along  with  a  more  detailed 
characterization  of  the  neutral-gain  utilization  for  variation  in  y0,  ®m,  and  ion  species  composition. 


C.  Decoupling  Voltage  and  Current  Utilization  from  Energy  Efficiency 

Energy  efficiency  is  factored  in  Eq.  (13)  as  the  product  of  voltage  utilization  and  current  utilization.  It  is 
convenient  to  introduce  %,  which  is  defined  in  Eq.  (14)  as  the  ratio  of  output  moles  of  charge  to  input  moles  of 
propellant.  Although  %  cancels  in  the  energy  efficiency  term,  it  is  shown  in  Equations  (17)  and  (18)  to  reveal  the 
physical  nature  of  the  voltage  and  current  utilization  inherent  in  the  thrust  efficiency. 
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Equation  (15)  transforms  the  first  term  in  brackets  of  Eq.  (13)  into  an  expression  that  contains  explicit 
utilization  of  acceleration  potential  for  each  ion  species,  8V/Vd.  An  ExB  probe  measurement  of  ion  mass  flow 
fractions  and  estimation  of  ion  acceleration  potentials  would  enable  a  voltage  utilization  to  be  calculated  for  each  ion 
species.  The  ratio  of  average  particle  specific  kinetic  energy  to  average  ion  specific  kinetic  energy  in  Eq.  (15)  is 
shown  to  be  approximately  equal  to  <Dm  in  Eq.  (16).  The  bracketed  quantity  in  Eq.  (16)  is  very  close  to  unity  under 
all  reasonable  conditions  of  Hall  thruster  operation  due  to  the  low  value  of  y0  in  the  numerator,  since  y02  (1- 
<I>m)<lxl04  for  experimental  measurements  of  yo<0.05.46'47  Thus,  the  approximation,  <v2>,„/<vi2>,„  =  <f>m  is  accurate 
to  better  than  one  in  ten  thousand. 
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To  simplify  the  analysis,  all  ions  are  considered  to  be  created  in  the  same  zone  whose  length  is  small 
compared  to  the  acceleration  length,  such  that  Va=8V;«constant.  Previous  investigations  by  Kim4*  and  King49  found 
that  species  dependent  energy  to  charge  ratios  varied  by  tens  of  volts  using  different  types  of  energy  analyzers.  This 
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variation  in  energy  to  charge  ratio  is  less  than  10%  of  the  typical  discharge  voltage.  According  to  Hofer50,  the 
approximation  will  have  a  negligible  effect  on  accuracy  since  the  plume  is  predominantly  composed  of  singly 
ionized  xenon.  The  velocity  of  neutrals  and  each  ion  species  is  approximated  with  a  delta  function  distribution  of 
velocities,  such  that  the  ion  species  kinetic  energies  are  proportional  to  their  charge  and  the  ion  velocity  ratio 
magnitudes  are  lyyl=(Zy)1/2.  To  the  author’s  knowledge,  this  approximation  is  consistent  with  all  other  performance 
architectures.  Therefore,  the  most  probable  ion  acceleration  potential,  Va,  as  measured  with  a  retarding  potential 
analyzer  (RPA)  enables  the  first  term  in  brackets  of  Eq.  (13)  to  be  expressed  in  terms  of  the  average  voltage 
utilization  efficiency  and  x  in  Eq-  ( 18).  The  average  voltage  utilization  in  Eq.  (18)  compares  the  most  probable 
potential  gained  by  ions  with  the  applied  anode  potential,  and  is  unity  for  ionization  at  the  anode  face  where  ions  are 
accelerated  through  the  entire  anode  potential. 
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The  second  bracketed  term  in  Eq.  (13)  is  transformed  in  Eq.  (19)  into  an  expression  containing  the  current 
utilization  efficiency  and  x-  Current  utilization  efficiency  is  the  fraction  of  cathode  electron  flow  that  electrically 
neutralizes  the  accelerated  positive  ions  in  the  plume,  and  is  calculated  as  the  ratio  of  ion  beam  current  to  discharge 
current. 
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Electrons  recycled  to  the  anode  ionize  neutral  propellant  and  sustain  dissipative  plasma  processes.  The 
minimum  power  required  to  produce  a  given  ion  species  composition  is  expressed  in  Eq.  (21).  Based  on  the 
minimum  power  requirement,  the  minimum  recycled  electron  current  in  a  Hall  thmster  discharge  may  be  determined 
as  a  function  of  the  applied  voltage  between  the  anode  and  cathode.  The  minimum  recycled  electron  current  fraction 
in  Eq.  (22)  shows  that  increased  discharge  voltage,  and  hence  increased  recycled  electron  energy,  reduces  the 
minimum  recycled  electron  current.  Theoretically,  a  stream  of  single  electrons  may  deposit  the  minimum  power  to 
produce  a  cascade  of  ionization.  Power  from  recycled  electrons  in  excess  of  Pmin  is  lost  by  Joule  heating  processes 
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to  the  plasma  and  channel  walls.  When  combined  with  the  total  power  loss  estimates  from  energy  efficiency,  this 
places  an  upper  limit  on  Joule  heating  losses  to  the  walls  amounting  to  (l-r|E-rmin)Pd. 
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D.  Evaluating  Plume  Divergence  Losses  in  Beam  Efficiency 

The  momentum-weighted  average  plume  divergence  is  defined  in  Eq.  (23)  as  the  ratio  of  the  measured 
thmst  component  directed  along  the  thruster  centerline  to  the  theoretical  thrust  achieved  when  all  ions  are  traveling 
parallel  to  thruster  centerline.  Momentum  losses  associated  with  plume  divergence  may  be  calculated  with 
knowledge  of  the  input  mass  flow,  measured  thrust,  and  the  mass-weighted  average  velocity. 
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The  value  <cos(9)>2  has  been  used  in  previous  analyses  to  describe  plume  focusing26,  but  there  has  not 
been  a  consistent  method  to  calculate  the  effect  of  plume  momentum  divergence  on  thrust.  This  is  primarily  due  to 
the  difficulty  of  measuring  particle  velocity  throughout  the  plume.  Charge  divergence  in  the  plume  is  indicative  of 
off-axis  velocity  losses  in  thrust,  and  is  a  useful  alternative  for  experimental  characterization  of  performance  losses 
due  to  plume  divergence.  The  momentum-weighted  average  divergence  is  approximated  as  the  charge -weighted 
average  divergence  in  Eq.  (24)  for  an  axisymmetric  plume,  which  enables  calculation  of  off-axis  cosine  losses  in  Eq. 
(25)  using  the  ratio  of  the  axial  component  of  beam  current  to  total  beam  current  as  measured  by  a  Faraday  probe. 
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Off-axis  cosine  losses  integrated  in  the  numerator  of  Eq.  (24)  quantify  the  axial  component  of  beam  current 
that  generates  thrust.  This  formulation  has  been  used  in  past  analyses24'51,  and  creates  a  method  where  the  plume 
divergence  vector  loss  is  evaluated  in  a  scalar  form.  The  calculation  of  is  complicated  by  the  presence  of  charge 
exchange  in  the  plume,  which  increases  measured  beam  current  at  large  angles  from  centerline  and  artificially 
increases  divergence  losses. 

An  effective  plume  divergence  angle  may  be  calculated  as  shown  in  Eq.  (26).  This  angle  is  significantly 
less  than  the  typically  reported  95%  divergence  half-angle. 
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The  ratio  m  (0)v(0)/J(9)  introduced  in  Eq.  (24)  reveals  the  difference  between  momentum-weighted 
divergence  and  charge-weighted  divergence.  This  ratio  is  evaluated  in  Appendix  C  for  variations  in  mass  utilization 
and  ion  species  population.  In  Eq.  (27),  the  ratio  is  evaluated  at  angular  position  0  and  reduced  to  a  function  of  Q 
and  <t>m. 
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Spatial  variation  of  the  average  ion  charge  and  mass  utilization  in  the  plume  is  the  primary  source  of 
discrepancy  between  momentum-weighted  divergence  and  charge -weighted  divergence.  Significant  variations  in  Q 
and  <Dm  within  the  angular  region  of  highest  beam  current  will  have  the  greatest  effect  on  plume  divergence 
calculations.  However,  these  variations  are  negligible  in  typical  Hall  thruster  and  ion  engine  plumes.  Variations  in 
propellant  efficiency  and  average  ion  acceleration  potential  in  the  plume  may  also  result  in  disparities,  but  are 
secondary  effects  and  not  expected  to  fluctuate  at  the  angular  region  of  peak  beam  current. 


E.  Evaluating  Exchange  Parameters  and  Thruster  Performance 

Two  experimental  parameters,  denoted  E|  and  E2,  may  be  written  using  the  preceding  utilization 
efficiencies  to  isolate  changes  in  performance  due  to  processes  related  to  thrust  or  discharge  current.  The 
parameters  are  formed  in  Eq.  (28)  such  that  the  product  of  Et  and  E2  is  equal  to  anode  thruster  efficiency.  Equations 
(29)  and  (30)  show  how  these  quantities  are  calculated  based  solely  on  telemetry  and  thrust  measurements,  and 
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provide  insight  about  the  relative  magnitudes  of  the  individual  utilization  efficiencies  in  the  absence  of  plume 
measurements. 
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For  fixed  thruster  telemetry  inputs  Vd  andrh,  the  dimensionless  experimental  parameter  Ei  separates 
changes  in  anode  efficiency  directly  to  variations  in  thrust  and  E2  isolates  changes  due  to  variations  in  discharge 
current.  The  quantity  E,~T2,  and  is  a  function  of  propellant  efficiency,  beam  efficiency,  voltage  utilization,  and  x- 
Experimental  parameter  Ei  relates  the  applied  acceleration  potential  to  dispersion  and  divergence  of  the  jet.  The 
quantity  E2~Id  *,  and  is  a  function  of  current  utilization  and  x-  Experimental  parameter  E2  relates  the  input  flow  of 
mass  to  the  total  output  flow  of  charge.  The  inverse  of  E2  was  used  in  the  Soviet  literature52  as  early  as  1978  and 
termed  the  exchange  parameter.53'54  The  naming  convention  used  here  describes  the  exchange  of  applied  input 
parameters  to  operational  thruster  properties.  Thus,  E|  is  termed  the  voltage  exchange  parameter  and  E2  is  termed 
the  mass  exchange  parameter.  While  ionization  and  acceleration  processes  are  closely  coupled,  the  form  of  the 
experimental  parameters  indicates  that  propellant  efficiency,  beam  efficiency,  and  voltage  utilization  are  principal  in 
the  formation  of  directed  thrust.  In  the  absence  of  diagnostics  for  the  determination  of  plasma  properties,  these 
experimental  parameter  groups  allow  limits  to  be  placed  on  acceptable  values  for  the  average  charge,  mass 
utilization,  and  plume  divergence. 55,56 

The  loss  mechanisms  that  effect  T/P  and  Isp  are  shown  in  Eq.  (31)  and  (32)  in  terms  of  the  experimental 
parameters  E,  and  E2  and  in  terms  of  the  utilization  efficiencies.  These  formulations  indicate  low  current 
utilization  and  large  x  will  decrease  T/P.  Conversely,  specific  impulse  is  not  directly  affected  by  current  utilization 
and  increases  for  large  x-  It  is  interesting  to  note  the  mass  utilization  in  propellant  efficiency  and  x  cancels  in  the 
formulation  of  T/P,  which  suggests  the  increased  T/P  at  low  discharge  voltage  is  not  directly  affected  by  ionization 
losses.57 


14 


The  T/P  formulation  is  similar  to  the  one  dimensional  analysis  used  to  estimate  thrust  of  the  SERT-II  ion 
engines,  which  is  equivalent  to  Eq.  (31)  when  the  term  (®P  ¥3  /  %)  is  unity.58-59-60-61-62’63  An  analogous  expression 
was  used  to  estimate  thrust  of  the  NSTAR  ion  engines  on-board  NASA’s  Deep  Space  1  by  including  factors  to 
estimate  the  effect  of  multiply  charged  ions  and  beam  cosine  losses.64  These  factors  were  similar  those  in  the 
formulation  of  thrust  by  Masek  et  al.  and  are  related  to  the  ratio  (®P  Tb  /  x)1/2  for  a  bimodal  ion  population.  17  The 
formulations  of  T/P  and  Isp  in  Eq.  (31)  and  (32)  will  be  used  in  addition  to  the  anode  efficiency  to  evaluate  the 
fidelity  of  plume  measurements  in  evaluating  thruster  performance. 


IV.  Application  to  Experimental  Results 

A.  Experimental  Apparatus 

The  efficacy  of  the  performance  architecture  is  studied  using  experimental  results  from  a  nominal  6  kW 
Hall  thruster  in  Chamber  3  at  the  Air  Force  Research  Laboratory.65  Operation  of  the  laboratory  model  thruster  with 
a  center-mount  cathode  ranges  from  105  V  to  300  V  with  10  mg/s  Xe  anode  flow.  Cathode  flow  fraction  (CFF  = 
cathode  flow  rate  /  anode  flow  rate)  is  increased  for  low  discharge  voltage  operation  to  maintain  stability  and 
maximize  total  efficiency. 

The  study  consists  of  thrust  measurements  and  examination  of  the  ion  voltage  distribution,  plasma 
potential,  ion  species  composition,  and  ion  current  density  in  the  far-field  plume.  Facility  effects  on  thrust  and 
current  density  profiles  are  minimized  by  extrapolating  measurements  at  several  background  pressures  to  zero 
pressure.  This  technique  is  commonly  used  to  approximate  thrust  in  the  space  environment  from  ground 
measurements,  and  was  recently  adapted  to  estimate  current  density  profiles  in  the  plume.66  Faraday  probe  sweeps 
were  executed  in  a  1  m  hemispherical  arc  from  0°  to  180°,  and  RPA  measurements  were  taken  at  several  angular 
positions  in  a  1  m  hemispherical  arc.  The  RPA  most  probable  ion  pass  voltage  was  corrected  with  Langmuir  probe 
plasma  potential  measurements  less  than  3  cm  from  the  probe,  which  resulted  in  negligible  variation  in  the  most 
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probable  ion  acceleration  potential  over  a  large  angular  span.  The  ion  species  composition  was  evaluated  from  1  m 


to  1.3  m  downstream  of  the  exit  plane  on  channel  centerline  and  current  fractions  were  corrected  for  the  loss  of 
beam  ions  due  to  CEX  collisions  using  the  procedure  outlined  by  Shastry.67  Conclusions  from  a  recent  study  of 
angularly  resolved  ExB  probe  spectra  using  this  thruster  model  asserted  that  a  single  point  channel  centerline 
measurement  was  representative  of  species  fractions  the  plume,  and  deviated  from  plume  averaged  results  by  less 
than  1.5%  over  a  wide  range  of  thruster  operating  conditions.68  ExB  and  RPA  results  in  the  next  section  are  not 
plume  averaged,  and  variations  are  accounted  for  in  the  specified  uncertainty.  Further  details  concerning  the 
facility,  diagnostics,  and  measurement  uncertainty  are  described  elsewhere.65 


B.  Experimental  Results  and  Analysis 

Experimental  results  and  measurement  uncertainty  are  listed  in  Table  1  for  constant  10  mg/s  anode  flow 
operation  from  105  V  to  300  V  discharge.  This  set  of  thrust  and  plume  measurements  enable  full  characterization  of 
the  thruster  utilization  efficiencies  and  performance  parameters. 


Table  1:  Thrust,  Telemetry,  and  Far-field  Plume  Measurements  of  a  6  kW  Laboratory  Hall  Thruster  65 


105  V,  10  mg/s 

16%  CFF 

120  V,  10  mg/s 

12%  CFF 

150  V,  10  mg/s 

7%  CFF 

300  V,  10  mg/s 

7%  CFF 

vd  [V] 

105  ±  0.05% 

120  ±  0.05% 

150  ±0.05% 

300  ±  0.05% 

Id  [A] 

9.8  ±  0.2% 

9.5  ±  0.2% 

9.1  ±0.2% 

8.8  ±  0.2% 

rh  [mg/s] 

10.0  ±  1% 

10.0  ±  1% 

10.0  ±  1% 

10.0  ±  1% 

T  [mN] 

87  ±  1% 

92  ±  1% 

112  ±1% 

184  ±  1% 

0.91  ±0.04 

0.90  ±  0.04 

0.86  ±  0.04 

0.82  ±  0.04 

Cl2 

0.09  ±  0.02 

0.10  ±0.02 

0.10  ±0.02 

0.12  ±0.02 

o3 

0.0  ±0.01 

0.0  ±0.01 

0.04  ±0.01 

0.06  ±  0.01 

Va  [V] 

81  ±8 

92  ±9 

120  ±9 

266  ±  10 

iBeam  [A] 

7.1  ±5% 

6.9  ±  5% 

7.1  ±5% 

7.4  ±  5% 

lAxial  [A] 

5.9  ±  5% 

5.8  ±  5% 

6.5  ±  5% 

7.2  ±  5% 

The  uncertainty  in  thmst  and  plume  measurements  were  propagated  into  calculations  of  plasma  properties 
and  utilization  efficiencies.  These  values  are  calculated  using  standard  error  propagation  techniques  (%  Error 
=  ^(tycErrort)2  +  (%Error,)2  +... )  for  independent  measurements.  Uncertainty  in  calculations  with  correlated  terms  is 
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conservatively  estimated.  Correlated  terms  include  the  axial  and  total  ion  beam  current  from  a  Faraday  probe  and 
the  ion  current  fractions  from  the  ExB  probe. 

Specific  impulse,  T/P,  and  anode  efficiency  based  on  plume  measurements  are  compared  to  the  values 
determined  by  thrust  measurements  in  Fig.  1 .  Agreement  between  the  performance  parameters  is  within  the  error 
bars  for  all  operating  conditions,  with  the  greatest  difference  occurring  at  105  V  and  120  V.  The  consistent 
agreement  between  thrust  measurements  and  plume  measurements  indicates  the  efficiency  architecture  is  accurately 
capturing  the  thmster  performance  characteristics.  Comparison  of  the  interrelated  T/P  and  Isp  terms  at  each 
operating  condition  support  the  measurement  validity  and  corroborate  the  formulation  of  loss  mechanisms  in  Eq. 
(31)  and  (32). 


Specific  Impulse  [  s  ] 

Fig.  1  Thrust  to  power  ratio  as  a  function  of  specific  impulse  with  lines  of  constant  anode  efficiency  for  a  6  kW  laboratory 
Hall  thruster  ranging  from  105  V  to  300  V  at  10  mg/s  operation. 


Anode  thmster  efficiency  is  separated  into  utilization  efficiencies  in  Fig.  2.  Reduced  anode  efficiency  at 
low  discharge  voltage  is  primarily  due  to  increased  beam  divergence  and  increased  energy  losses.  Beam  efficiency 
decreases  from  0.94  at  300  V  to  0.69  at  105  V,  which  corresponds  to  an  increase  in  divergence  half-angle  from  14° 
to  34°.  The  decline  in  energy  efficiency  with  discharge  voltage  is  the  dominant  loss  mechanism  during  low  voltage 
operation.  Propellant  efficiency  varies  by  less  than  2%  over  the  range  of  operation,  and  shows  a  slight  increase  at 
the  lowest  discharge  voltages.  This  effect  is  studied  by  decomposing  propellant  efficiency  into  mass,  charge,  and 
neutral-gain  utilization  efficiencies  in  Fig.  3. 
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The  charge  utilization  decreases  at  higher  discharge  voltage  due  to  the  production  of  multiply  charged  ions, 
which  is  consistent  with  past  investigations.43'6*  Increased  mass  utilization  at  low  voltage  is  most  likely  caused  by 
additional  ionization  as  a  consequence  of  the  increased  CFF,  and  is  responsible  for  the  increased  propellant 
efficiency  at  low  voltage.  Variations  in  neutral-gain  utilization  due  to  neutral  speed  are  quantified  in  Appendix  B. 
The  neutral  speed  is  conservatively  estimated  at  200  m/s  based  on  LIF  studies  of  this  thruster  model  at  the 
University  of  Michigan.47 


Discharge  Voltage  [V] 

Fig.  2  Anode  efficiency  and  utilization  efficiencies  as  a  function  of  discharge  voltage  for  a  6  kW  laboratory  Hall  thruster 
ranging  from  105  V  to  300  V  at  10  mg/s  operation. 


Discharge  Voltage  [VI 

Fig.  3  Propellant  efficiency,  mass  utilization,  charge  utilization,  and  neutral-gain  utilization  as  a  function  of  discharge  voltage 
for  a  6  kW  laboratory  Hall  thruster  ranging  from  105  V  to  300  V  at  10  mg/s  operation. 
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Losses  in  energy  efficiency  are  divided  into  current  utilization  and  voltage  utilization  in  Fig.  4.  At  low 
voltage,  the  ionization  cost  becomes  a  larger  faction  of  the  applied  anode  potential.  The  decrease  in  voltage 
utilization  is  approximately  equal  to  the  fractional  increase  of  the  xenon  ionization  potential  relative  to  the  thruster 
discharge  voltage  (12  eV  /  Vd).  The  fraction  of  electron  current  to  the  anode  is  the  dominant  loss  mechanism  in 
energy  efficiency,  and  is  magnified  during  low  voltage  operation. 


Fig.  4 


Discharge  Voltage  [V] 

Energy  efficiency,  voltage  utilization,  and  current  utilization  as  a  function  of  discharge  voltage  for  a  6  kW  laboratory 
Hall  thruster  ranging  from  105  V  to  300  V  at  10  mg/s  operation. 


In  Fig.  5,  the  fraction  of  recycled  electron  current  is  compared  to  the  minimum  electron  current  to  the 
anode  necessary  to  sustain  ionization.  As  thruster  power  is  decreased  for  constant  mass  flow  operation,  the 
measured  anode  electron  current  fraction  increased.  The  increase  in  r  was  larger  than  the  increase  in  rmln,  and 
highlights  a  significant  loss  mechanism  associated  with  low  voltage  Flail  thruster  operation.  The  decrease  in  rmin 
with  increased  discharge  power  is  expected  for  higher  discharge  voltage  operation  from  Eq.  (22). 

Total  power  losses  as  a  function  of  discharge  power  are  estimated  in  Fig.  6.  The  minimum  power  loss  to 
sustain  ionization,  Prnln  -120  W,  is  nearly  constant  from  1000  W  to  2600  W  discharge  power.  Total  Joule  heating 
from  cathode  electrons,  Pdrmin,  is  relatively  constant  below  150  V  discharge,  and  increases  by  50%  to  300  W  during 
300  V  operation.  This  increased  Joule  heating  from  recycled  electrons  with  increased  discharge  power  indicates 
excessive  losses  to  the  channel  walls  and  plasma  excitation  processes.  The  total  power  loss  in  the  thruster  is 
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estimated  using  energy  efficiency.  It  is  important  to  note  these  power  losses  do  not  account  for  the  spread  in  most 
probable  ion  potential  or  variations  in  the  plume,  which  introduce  significant  error  in  the  estimates.  Further  isolation 
of  energy  loss  mechanisms  are  unfeasible  with  plume  measurements,  and  require  numerical  simulations  to  determine 
the  effects  of  wall  collisions,  radiation,  excitation  and  relaxation  processes. 


Discharge  Power  [  W  ] 

Fig.  5  Fractional  anode  electron  current  and  minimum  anode  electron  current  as  a  function  of  discharge  power  for  a  6  kW 
laboratory  Hall  thruster  ranging  from  105  V  to  300  V  at  10  mg/s  operation. 


Fig.  6 


Discharge  Power  [  W  ] 

Power  losses  as  a  function  of  discharge  power  for  a  6  kW  laboratory  Hall  thruster  ranging  from  105  V  to  300  V  at  10 
mg/s  operation. 


20 


V.  Summary 


The  separation  of  scalar  thrust  into  mass-weighted  and  momentum-weighted  terms  enabled  factorization  of 
anode  efficiency  into  the  product  of  (1)  energy  efficiency,  (2)  propellant  efficiency,  and  (3)  beam  efficiency.  The 
analysis  decomposed  anode  efficiency  from  first  principles,  and  formulated  the  relationship  such  that  losses 
associated  with  energy  conversion  are  analyzed  separately  from  losses  associated  with  dispersion  of  the  jet  VDF  and 
beam  divergence.  The  effects  of  multiply  charged  ions  are  included,  and  the  neutral-gain  utilization  is  introduced  to 
account  for  the  neutral  speed.  The  momentum-weighted  divergence  loss  term  (<cos(0)>2mv)  contains  all  jet  vector 
losses,  and  differences  from  the  current-weighted  divergence  was  characterized.  The  proposed  analysis  is  compared 
to  past  performance  methodologies,  with  similarities  to  previous  architectures  formulated  by  Masek  et  al.  in  the  ion 
engine  community17,  Belan  et  al.  from  the  Kharkov  Aviation  Institute24,  and  a  contemporary  methodology  by 
Hofer.43'50 

A  case  study  was  presented  with  experimental  performance  and  far-field  plume  measurements  for  a  6  kW 
laboratory  Hall  thruster.  Comparison  of  calculated  anode  efficiency,  T/P,  and  Isp  for  constant  anode  flow  operation 
from  105  V  to  300  V  showed  good  agreement  between  the  analytical  model  and  plume  diagnostics.  The 
performance  methodology  highlights  low  current  utilization  and  beam  divergence  as  the  primary  loss  mechanisms 
during  low  voltage  operation.  Current  utilization  can  never  be  unity  due  to  a  finite  flow  of  electrons  that  are 
recycled  to  the  anode  to  sustain  ionization  processes.  Power  losses  and  Joule  heating  processes  are  estimated  for 
each  operating  condition.  The  minimum  power  loss  to  sustain  ionization  showed  minimal  variation  with  increasing 
discharge  power  for  constant  anode  mass  flow  operation.  For  increasing  voltage  utilization  with  discharge  voltage, 
the  increasing  power  losses  during  high  voltage  operation  are  the  result  of  additional  wall  collisions,  radiation, 
excitation  and  relaxation  processes  within  the  discharge  channel. 
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APPENDIX  A:  Definitions  of  Mass-weighted  and  Momentum-weighted  Quantities 

In  Eq.  (2),  the  steady  state  equation  for  axisymmetric  thrust  directed  along  thruster  centerline  was  factored 
into  the  product  of  mass  flow,  mass-weighted  average  velocity,  and  momentum-weighted  average  divergence  using 
the  definitions  shown  in  Eq.  (A-l),  (A-2),  (A-3),  and  (A-4). 

Hemispherical  integration  of  mass  flow  rate  throughout  the  plume  at  constant  radius  in  Eq.  (A-l)  is  equal  to 
the  total  mass  flow  rate  supplied  to  the  anode.  Propellant  flow  exiting  the  cathode  has  a  negligible  effect  on  thmst 
since  it  is  composed  primarily  of  neutrals.  Thus,  cathode  flow  is  considered  a  parasitic  loss  and  is  not  included  in 
this  definition. 

31/2 

hi  =  2 jrR2  j* iii(9) sin(0) d9  (A-l) 

0 

The  average  velocity  vector,  v  ,  is  defined  in  Eq.  (A-2)  as  the  integral  of  the  VDF  over  all  velocity  space  at 
for  an  axisymmetric  plume  at  angular  position  9  and  radius  R.  The  radial  propellant  velocity  component  in  spherical 
coordinates,  v  ,  is  formulated  in  Eq.  (A-3)  and  radial  ion  velocity  component,  vj ,  is  formulated  in  Eq.  (A-4)  using 
ion  flow  fractions.  This  is  the  dominant  velocity  component  for  ions  originating  from  the  thruster  discharge,  and 
other  velocity  components  in  the  plume  are  negligible.  Squared  ion  and  propellant  radial  velocity  terms  are  found 
using  analogous  expressions. 


v=JjJvf(v)dv/JJJf(v)dv 


(A-2) 


V 


mj 
4-i  m 


j= 0 


(A-3) 


n =y— v;  =yfj 


(A-4) 


The  mass-weighted  average  radial  velocity  component  for  an  axisymmetric  plume  is  defined  in  Eq.  (A-5) 
and  simplified  using  Eq.  (A-l).  The  momentum-weighted  average  cosine  is  defined  in  Eq.  (A-6)  and  simplified 
using  Eq.  (A-5).  The  definition  of  thmst  in  Eq.  (2)  is  naturally  produced  through  the  relation  of  mass-weighted  and 
momentum-weighted  terms  in  Eq.  (A-6).  This  formulation  of  thrust  sets  the  foundation  for  the  Hall  thmster 
efficiency  architecture  outlined  the  paper. 
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All  quantities  are  assumed  steady  state  and  evaluated  at  constant  radius  from  the  exit  plane,  such  that  the 
quantities  to  be  integrated  have  units  of  sr'1.  For  hemispherical  integration  about  a  point  source,  the  measurement 
distance  must  be  multiple  thruster  diameters  downstream  of  the  exit.  The  definitions  in  Eq.  (A-l)  through  (A-6)  do 
not  account  for  dispersion  of  the  jet  due  to  scattering  and  charge  exchange  ions  in  the  plume.  These  processes  will 
have  a  different  effect  depending  on  the  diagnostic,  facility,  background  pressure,  and  distance  from  the  thruster.  To 
this  point,  an  axisymmetic  plume  with  radial  particle  velocity  and  negligible  facility  effects  are  the  only 
approximations  required. 
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APPENDIX  B:  Formulation  of  Neutral-Gain  Utilization  in  Propellant  Efficiency 

Propellant  efficiency  characterizes  dispersion  of  the  VDF  and  is  a  measure  of  the  jet  momentum  relative  to 
the  jet  kinetic  energy.  In  Eq.  (B-l)  and  (B-2),  the  second  term  in  brackets  in  Eq.  (7)  is  decomposed  into  the  product 
of  mass  utilization  and  neutral-gain  utilization. 
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The  neutral-gain  utilization  is  defined  in  Eq.  (B-3),  and  is  always  greater  than  unity.  The  approximation  in 
Eq.  (12)  is  formulated  by  neglecting  second  order  terms  with  y0,  and  shows  <bN_G  is  primarily  a  function  of  the 
neutral  speed  and  mass  utilization.  Effects  that  increase  the  neutral-gain  utilization  result  in  detrimental  losses  to 
other  utilization  efficiencies.  Ideal  thruster  operation  would  correspond  to  unity  neutral-gain  for  100%  ionization  or 
zero  neutral  speed  for  infinite  neutral  residence  time  in  the  channel. 
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The  small  gain  in  efficiency  due  to  the  neutral  speed  is  generally  neglected.  Neutral-gain  utilization  is 
plotted  in  Fig.  (B-l)  for  Q=l,  and  is  shown  to  increase  with  reduced  mass  utilization  and  large  neutral  speed. 
Thrusters  with  a  high  neutral  thermal  speed  or  large  axial  injection  speed  may  exhibit  neutral-gain  utilization  of  1.01 
to  1.02.  In  Fig.  (B-2),  the  neutral-gain  utilization  is  shown  for  the  ion  species  composition  of  the  6  kW  laboratory 
thruster  studied  in  Section  IV.  Mass  utilization  of  Hall  thrusters  typically  ranges  from  0.80  to  0.90  for  nominal 
operation.  For  estimates  of  yo®0.02,47  the  neutral-gain  utilization  is  approximately  1.005  when  mass  utilization  is 
0.90.  The  neutral-gain  is  expected  to  be  larger  for  low  discharge  voltage,  high  power  operation  where  the  value  of 
y0  is  increased  and  ionization  is  reduced. 
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Mass  Utilization,  <t>m 

Fig.  B-l  Neutral-gain  utilization  as  a  function  of  mass  utilization  with  variations  in  reduced  neutral  speed  from  y()=0.01  to 
yo=0.O6  for  Q=l. 
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Fig.  B-2  Neutral-gain  utilization  for  the  ion  species  composition  of  the  6  kW  laboratory  Flail  thruster  during  300  V,  10  mg/s 
operation  as  a  function  of  mass  utilization  for  y0=0.01,  0.03,  and  0.06.  Neutral-gain  of  the  experimental  ion 
composition  is  bounded  by  lines  of  constant  Q=1  and  Q=3,  which  are  limiting  cases  for  a  trimodal  ion  population. 
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APPENDIX  C:  Comparison  of  Momentum-weighted  and  Charge- weighted  Plume  Divergence 


The  momentum- weighted  average  divergence  is  approximated  as  the  charge -weighted  average  divergence 
in  Eq.  (24),  which  enables  calculation  of  off-axis  cosine  losses  using  the  ratio  of  the  axial  component  of  beam 
current  to  total  beam  current  as  measured  by  a  Faraday  probe. 

Analysis  of  the  ratio  m  (9)v(0)/J(9)  introduced  in  Eq.  (24)  characterizes  the  properties  that  cause 
differences  between  momentum-weighted  divergence  and  charge-weighted  divergence.  The  ratio  is  evaluated  at 
angular  position  9  In  Eq.  (C-l).  Each  particle  species  is  approximated  with  a  delta  function  velocity  distribution  and 
velocity  ratios  are  calculated  for  the  idealized  case  of  ion  creation  at  the  same  location.  The  ratio  m  (9)v(9)/J(0)  is 
characterized  using  the  dimensionless  quantity  y,  which  describes  differences  associated  with  ionization  fraction  and 
ion  species  population.  The  term  y  is  shown  as  a  function  of  <tim,  <tiP,  and  y0  in  Eq.  (C-2). 
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The  term  inside  the  brackets  of  Eq.  (C-2)  is  near  unity  for  yo=0.02,47  and  variations  in  the  neutral  speed 
ratio  have  negligible  effect  on  y.  In  Figure  C-l,  y  is  shown  for  mass  utilization  ranging  from  0.6  to  1.0  and  yn=0.02. 
The  ratio  is  bounded  for  a  trimodal  ion  population,  and  is  shown  to  be  primarily  a  function  of  mass  utilization  and 
average  charge. 

Variations  in  the  magnitude  of  <I>m  and  Q  alter  y,  but  the  angular  location  of  these  variations  is  the  dominant 
factor  causing  differences  between  momentum-weighted  divergence  and  charge -weighted  divergence.  This  effect  is 
illustrated  in  Figure  C-2,  which  shows  a  representative  distribution  of  current  density  and  beam  current  as  a  function 
of  angular  position  for  the  6  kW  laboratory  Hall  thruster  from  Section  IV  operating  at  300  V.  The  location  of  peak 
beam  current  in  the  plume  is  approximately  8  degrees  from  thruster  centerline,  and  variations  in  <£>m  and  Q  within  the 
angular  range  of  full  width  half  maximum  (FWHM)  will  have  the  greatest  effect  on  divergence.  For  the  beam 
current  in  Figure  C-2,  this  range  is  approximately  ±3  to  ±21  degrees.  Near  and  far-field  plume  measurements  of  the 
the  6  kW  Hall  thruster68,  the  SPT-10048'69  and  the  BHT-20070  Hall  thruster  found  that  the  ion  species  populations 
did  not  significantly  change  within  +15  degrees  from  thruster  centerline,  and  the  fraction  of  Xe+2  increased  sharply 
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by  approximately  5%  at  ±20  degrees.  These  changes  in  ion  species  composition  have  a  negligible  effect  on  the 
value  of  <cos(9)>mv,  thus  enabling  the  approximation  of  equivalence  between  momentum-weighted  divergence  and 
charge-weighted  plume  divergence  as  measured  with  a  Faraday  probe. 


Fig.  C-l  Variation  in  y  due  to  mass  utilization  and  ion  species  composition.  Mass  utilization  manifolds  are  shown  for  ®m  =  1.0, 
0.9,  0.8,  0.7,  and  0.6.  Lines  of  constant  f2  =0  and  f3  =0  bound  y  for  a  trimodal  ion  population,  and  are  compared  with 
the  approximation  y  ~  ®mQ'1/2. 
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Angular  Position  [Degrees] 


Fig.  C-2  Representative  distribution  of  beam  current  and  current  density  at  1  meter  radius  as  a  function  of  angular  position  in  the 
plume  of  a  nominal  6  kW  laboratory  Hall  thruster  operating  at  300  V,  10  mg/s  from  Section  IV. 
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